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We obtain the symmetry algebra of multi-matrix models in the planar 
large N c limit. We use this algebra to associate these matrix models with 
quantum spin chains. In particular, certain multi-matrix models are exactly 
solved by using known results of solvable spin chain systems. 



1 



Quantum systems whose degrees of freedom are matrices appear in several 
areas of mathematics and physics; for example, Yang- Mills theory H, ||, 
, string theory || , Q and M -theory {?j , || . Of particular interest is the 
limit as the dimension N c of the matrices goes to infinity. In this limit the 
dynamics is expected to simplify; for example, the quantum fluctuations of the 
invariants are of order 1/N C . The algebra of invariant observables becomes 
a Poisson algebra discovered in Ref. ||. For the general large N c limit, 
these Poisson brackets are very non-linear. The planar large N c limit is 
equivalent to a further approximation that replaces this Poisson algebra with 
a Lie algebra. In this paper we will describe this Lie algebra of observables 
of the matrix model in the planar limit, by a direct argument. 

As an illustration of the power of this new symmetry algebra, we will use 
it to solve some matrix models in the large N c limit. More precisely, we will 
map certain matrix models to quantum spin chains and use results from the 
theory of spin chains to solve them. This is reminiscent of the work || that 
connects some integrals over finite chains of matrices with classical integrable 
systems. From this point of view, our result is that certain path integrals over 
matrices can be mapped into quantum integrable systems. However we will 
mostly use the canonical formulation rather than the path integral formulation 
of these systems. 

We will study a class of matrix models whose degrees of freedom are a set 
of matrix-valued bosonic variables a^(i), aJ, M (i) satisfying the canonical com- 
mutation relations [a^(i),aP(j)] = [al tJ, (i),ajP(j)]=0 and [a^(i), <4- p (j)] = 
8(i,j)5£5P. Here, /i, v = 1, 2, . . . or N c . The positions of the indices 
indicate the transformation properties under U(N C ) : a 1 ^ —> g^g*" a p a etc. 
The degree of freedom labelled by the indices /U, u etc. will be called 'color' 
in analogy with quantum chromo dynamics (QCD). Indeed our matrix model 
can be thought of as a regularized version of pure QCD, with the variables 
a, a) representing gluons. The indices i = l,---,A describe the degrees of 
freedom (other than color) of the system. The Hamiltonian (along with all 
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other observables) will be required to be color invariant, i.e., invariant under 
the adjoint action of U(N C ) on a and a) . 

The path integral over matrix- valued functions of time, Ptf(j, t), Q£(j, t) 
with Lagrangian 

L(P,Q) = Y,P^t)j t Qim - H(P(t),Q(t)) 
j=i 

gives an equivalent theory, with the identifications = QH(j) + 

^Pv{j)i a u (J) = QuU) ~ iPv{j) j but the canonical formulation is more con- 
venient for our purposes. 

Define the vacuum state of the representation of these relations by a|0 >= 
. In the limit of large N c the color invariant states of the system are the 
'closed string' (or 'glueball') states such as 

= N- c/2 a^{k x )a^{k 2 ) ■ ■ ■ a^(k c )\0). 

Here strings of indices are denoted by capital letters. For example, K 
stands for k%, ■ ■ ■ , k c . The state is invariant under cyclic permutations; the 
equivalence class of permutations related to K by cyclic permutations is 
denoted by (K) . 

The operators that dominate the large iV c limit are 

gIj s N-^ b - 2 ^a^(h)a^(i2) ■ ■ ■ afriia) 

(Notice the reversal of order in the indices in the string J ; this definition 
serves to simplify some later equations.) All observables of a matrix model 
which survive in the large N c limit- the Hamiltonian of regularized QCD 
for example - are linear combinations of such operators. These states and 
operators were previously studied in Ref. 0, where an elegant application to 
large N c QCD is described . 

The factors of iV c have been chosen to obtain the 'planar' limit; it is 
so called because in perturbation theory, the Feynman diagrams that survive 
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can be drawn on a plane. There are other ways of taking the large N c limit, 
but the planar limit is the simplest. 

In the limit as N c — > oo these operators will map single closed string 
states to linear combinations of single closed string states ("glueballs") : 

This is the key simplification of the planar limit. (To higher orders in the 
1/N C expansion, there will be terms that correspond to splitting a glueball 
into several glueballs.) Here, SKs is equal to the number of different cyclic 
permutations of J such that each permuted sequence is identical to K. Also, 
in the second term we sum over all ways of splitting the sequence (K) into 
non-empty subsequences K\ and K% . A graphical representation of (|l]) 
is given in Fig. 1. 

The operators gj are like matrices except that they operate on the 
space of cyclically symmetric tensors. We will call them 'cyclix' operators. 
The product 9j9l °f two of the above operators is not a finite linear 
combination of the g 's themselves. But the commutator is indeed such a 
finite linear combination: finite linear combinations of the operators gj form 
a Lie algebra. (By finite linear cobinations we mean a sum over all sequences 
of indices I and J , of the form c i9j > such that only a finite number 
of the coefficients cj are non-zero.) The discovery of this Lie algebra is our 
main result. We will see that it has powerful consequences: for example we 
can solve some matrix models exactly using this newly discovered dynamical 
symmetry. 

Before we describe the commutation relations between two g l j 's , it 
is convenient to introduce another kind of operator fiQ on closed string 
states. The defining equation for these operators is fnl^ K ^ = ^(j)^^- 
These are thus the Weyl matrices in the basis of closed string states 

up to constant multiples. Rather than being independent operators, they 
are in fact just linear combinations of gj : f91 = 9j — J2k=i 9jk anc ^ 
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/(j) = 9 j ~ J2k=i 9k J- The two different ways of writing f91 imply that 
the operators gj are not linearly independent. 

Now we can state the commutation relations of our Lie algebra: 

'&&]=6Wl+ E E 

JiJ 2 =J K 1 K 2 =K 

+ E ML a + E + E W 

Ji.7 2 =J K 1 K 2 =J 

JlJ 2 = J 

K 1 K 2 = K 

+ E *3?«5f + E 

J1J2 Js=J 

JlJ 2 = J 
K X K 2 = K 

+ E « 1/X3 + E 

JlJ 2 = J 

K X K 2 = if 

+ E + E W$f 2) 

J1J2J3 = <^1<^2 = J 

K X K 2 = K K X K 2 K 3 = K 

J1J2J3 = J 
K X K 2 K 3 = K 

Although it appears complicated when written this way, these commutation 
relations have a rather natural graphical interpretation which we will describe 
in a longer paper |10|| . We will call the Lie algebra defined by these commu- 
tation relations the 'cyclix Lie algebra' or C_ M . 

The above defined ffy span an ideal of this algebra isomorphic to the 
inductive limit of linear algberas, gl^ . ( 5/00 can also be defined as the 
Lie algebra of matrices with only a finite number of non-zero entries.) 

We can quotient C_ M by this ideal to get another Lie algebra C_m > 
which is the essentially new object we have discovered. However it is only the 
extension C_ M that has a representation on the space of closed string states. 
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In the simplest special case of a matrix model with just one degree of 
freedom ( M = 1 ), the algebra C_y is just the algebra of (polynomial) 
vector fields on the circle. C_i is then the extension of this algebra by the 



algebra of finite-rank matrices [11|. Perhaps then C_m can be realized as 
the Lie algebra of vector fields on a non-commutative manifold. 

We will now show how some large N c matrix models can be solved 
by using this new symmetry algebra. Suppose the Hamiltonian of a matrix 
model is a linear combination H = J2u hjdj where hj = unless / and 
J have the same number of indices. ( This means that the 'gluon number' 
is a conserved quantity: regularized QCD is not of this type.) Such linear 
combinations form a subalgebra; let us call it C_ M . 

There is an isomorphism between multi-matrix models whose hamiltonians 
are in cf M and quantum spin chains. Now, there are some well-known 
examples of exactly solved quantum spin chains; they yield exactly solved 
matrix models. 

More explicitly, consider a spin chain with v sites: at any site a = 1 
,. . . , or v , there is a variable i a (called 'spin' for historical reasons) that 
can take the value 1,2, . . . , or A . We will impose the periodic boundary 
condition. A basis of states is given by \k\ ■ ■ ■ k v > . 

Define the operator 

X){a)\ki ■■■k v >= 5f \ki--- k a ^ik a+1 ■ ■ ■ k v > . 

This is just the Weyl matrix at site a . Let Aj(a) = X l Aa — v) if a > v. 
Now we can check that if / and J have the same length b < v , 

rM) ^E^! (a + l)---X;i(a + b-l) 

a=l 

satisfies the commutation relations of the algebra cf^j . If we also set 
r v(dj) = f° r b > v , we will have a representation r v of C_ M . The 
states of the periodic spin chain with zero total momentum correspond to 
cyclically symmetric tensors which are the states of the matrix model. 
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With each matrix model whose hamiltonian H = J2u hjdj is i n Cju 
, we can associate a quantum spin chain with the Hamiltonian 

F s P in = X ?M X h( a + 1) • • • Xf b {a + b - 1). 

IJ a=l 

Thus matrix models conserving the gluon number correspond to quantum 
spin systems with interactions involving neighborhoods of spins {a, a+1, a+ 
2,.. .,o + 6-l} . 

Let us look at some examples of solvable spin models and their associated 
matrix models. The simplest solvable quantum spin chain is perhaps the Ising 



model |1, [13]: 



= E ^» + A E t»t> + 1). 

a=l a=l 

Here r'^ ,y ' z are Pauli matrices at site a . Let the states 1 and 2 in the 
matrix model correspond to the spin-up and spin-down states in the Ising 
model. Using the fact that = X\ (a) - Xf (a) and + ir% = 1X\ (j), 
we get the corresponding element in : 

^Ising - 9l ~ 92 + A[9ll + 9l2 + 521 + 522 J- 

This is the large iV c limit of the matrix model with the Hamiltonian 

flSg*= tr [ot(l)a(l)-at(2)o(2)] + A tr [ a t(2)ot(2)a(l)a(l) 

+a t (2)a t (l)a(2)a(l) + a t (l)a t (2)a(l)a(2) + a t (l)a t (l)a(2)a(2)]. 



Our results, along with known results of the Ising spin chain [13] give the 
spectrum of this matrix model in the large N c limit: 

E(n p , u) = -2 £ (l + 2A cos h^r] + A 2 ] n p 
p=— f \ / 

where is any positive integer and n p = or 1. Also, we must impose the 

condition J2p=-v n pP = to get cyclically symmetric states. In particular 

we see that the value A = 1 is the critical value of the matrix model at 

which the spectum (in the planar limit) is that of a massless free fermion field 

on a lattice. 



It is interesting to ask whether the symmetries of the Ising spin chain can 



be understood within our formalism . Recall that [[12j|, [14] the solvability 
of the Ising model is due to the existence of an infinite number of conserved 
quantities. They form an infinite-dimensional Lie algebra, the Onsager al- 
gebra. This is the Lie algebra generated by iterating commutators of two 
operators Hq and V satisfying [Ho, [Ho, [Ho, V]]] = 16[Hq,V] and 
[V, [V, [V,H ]\] = 16[V,H }. For the Ising model, H = H = g{ - g\ and 
y = 5ii + 512 + 521 + 522 • Clearly, the Onsager algebra is a subalgebra of 
Cjyl . In particular, all conserved quantities of the Ising model are contained 
in our cyclix Lie algebra. It is not known whether this Ising matrix model is 
solvable for an arbitrary finite value of N c . 

To every solved spin chain there is thus a corresponding solved matrix 
model. Instead of a comprehensive list, we are just going to give a few illus- 
trative examples. 



The generalization of the Ising model with the Hamiltonian [15] 

a=l a=l 

also has the Onsager algebra as a dynamical symmetry. It corresponds to the 
element 

^matrix = g l _ g 2 + x[g 22 + (1 _ + (1 + 2iv)g X2 + ^ 

of the cyclix Lie algebra and hence to the exactly solvable matrix model 

Faf" x = tr [at(l)o(l)-ot(2)a(2)] + A tr [ t(2)at(2)a(l)a(l) 

+(1 - 2i«)a + (2)a + (l)a(2)a(l) + (1 + 2w)a t (l)a t (2)a(l)a(2) 
+a t (l)a t (l)a(2)a(2)]. 

The XYZ model fl]| with the Hamiltonian 

a=l a=l 

is a generalization of the Ising model in another direction. The corresponding 
element in the cyclix algebra is 
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mf x = all - oil - all + all - A[(i - + <&) + (1 + v)( g % + g %)]. 

A special case of this, the equivalence of a matrix model to the XXZ model, 
was found in |17[ |. 

The above correspondence between spin chains and matrix models is not 



restricted to the case «M = 2 . The chiral Potts model [18] has the Hamil- 
tonian 

a=l fc=l 

where afc,cifc are constants. Also, P a and Q a are generalized spin 
matrices at site a : Q = diag (1, ui, uj 2 , ■ ■ ■ ar^ -1 ) and P is defined 
by PQ = ujQP . Here, u> = e 2m ^ . This model is exactly solvable and 
corresponds to the element 

k=l j=l iij2=l 

where j\ + k should be replaced with ji + k — Ai£jx + k > A and j% — k should 
be replaced with j2 + A — k if J2 — k < in djljt^' : ' 2 ~ k of the cyclix algebra. 

The problem of finding the partition function of the spectrum of a hamil- 
tonian H is equivalent to evaluating the path integral over paths of period 
T : Tr e~ iHT = J D[P]D[Q]e* XCi tr mQU)-H{P(j),QmW _ wherg 

-fT(P, Q) is obtained by substituting a = Q + iP, = Q — iP into -ff as 
described previously. By applying this transcription to the above systems, we 
can obtain path integrals over matrices which can be evaluated exactly in the 
planar large N c limit. We wont give explicit expressions to keep the paper 
short. 

In addition to integrable matrix models associated with quantum spin 
chain models, we have also formulated models for QCD in terms of elements 
of the cyclix algebra [10|. We have also found the analog of the cyclix algebra 



suitable for studying open strings ('meson states') | |19| ; the supersymmetric 
extension has also been constructed [BfJ. The former is of interest in spin 
chains with open boundary conditions and QCD with quarks, and the latter 
in M -theory. 
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Iv 



J* K( J 




(a) (K) = (J) 





(b) Kj = J 

FIG. 1. The action of a gluonic operator on a single glueball state. The gluonic operator gj 
searches for a substring of K that agrees with J . If found, it replaces each such substring with 
/ ; otherwise, we get zero. Here, J* denotes the reverse of the sequence J . 
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